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ABSTRACT 

We study the thermal and chemical evolution during the Kelvin-Helmholtz 
phase of the birth of a neutron star, employing neutrino opacities that are 
consistently calculated with the underlying equation of state (EOS). Expressions 
for the diffusion coefficients appropriate for general relativistic neutrino transport 
in the equilibrium diffusion approximation are derived. The diffusion coefficients 
are evaluated using a field-theoretical finite temperature EOS that includes 
the possible presence of hyperons. The variation of the diffusion coefficients is 
studied as a function of EOS and compositional parameters. We present results 
from numerical simulations of protoneutron star cooling for internal stellar 
properties as well as emitted neutrino energies and luminosities. We discuss the 
influence of the initial stellar model, the total mass, the underlying EOS, and 
the addition of hyperons on the evolution of the protoneutron star and upon the 
expected signal in terrestrial detectors. 

We find that the differences in predicted luminosities and emitted neutrino 
energies do not depend much upon the details of the initial models or the 
underlying high-density EOS for early times (t < 10 s), provided that opacities 
are calculated consistently with the EOS. The same holds true for models which 
allow for the presence of hyperons, except when the initial mass is significantly 
larger than the maximum mass for cold, catalyzed matter. For times larger 
than about 10 seconds, and prior to the occurrence of neutrino transparency, 
the neutrino luminosities decay exponentially with a time constant that is 
sensitive to the high-density properties of matter. We also find the average 
emitted neutrino energy increases during the first 5 seconds of evolution, and 
then decreases nearly linearly with time. In general, increasing the protoneutron 
star mass increases the average energy and the luminosity of neutrinos, as well 
as the overall evolutionary time scale. The influence of hyperons or variations 
in the dense matter EOS is increasingly important at later times. Metastable 
stars, those with hyperons which are unstable to collapse upon deleptonization, 
have relatively long evolution times, which increase the nearer the mass is to the 
maximum mass supported by a cold, deleptonized star. 

Subject headings: protoneutron star, neutrinos, supernova, stellar evolution 
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1. INTRODUCTION 

A protoneutron (PNS) star forms in the aftermath of a successful supernova explosion 
as the stellar remnant becomes gravitationally decoupled from the expanding ejecta. The 
neutrinos radiated by the remnant are important to supernova energetics, and possibly 
are essential to supernova shock revival (Wilson 1985, Burrows, Hayes, & Fryxell 1995) 
They are also crucial to supernova nucleosynthesis, especially the r-process which may 
occur in the ejecta (Woosely et al. 1994). A careful study of the energies and time scales 
of neutrino emission in all flavors is also necessary for the interpretation of the data that 
will be recorded by terrestrial neutrino detectors in the event of a supernova occurring in 
our Galaxy or within several hundred kiloparsecs. In principle, observations of supernova 
neutrinos may not only allow an assessment of the remnant mass, but may also discriminate 
among competing types of equations of state (EOS) and internal neutron star compositions. 

The essential microphysical ingredients that govern the macrophysical evolution of the 
PNS in the so-called Kelvin-Helmholtz epoch, during which the remnant changes from a 
hot and lepton-rich PNS to a cold and deleptonized neutron star, are the (EOS) of dense 
matter and its associated neutrino opacity. Among the characteristics of matter that widely 
vary among EOS models are their relative compressibilities (important in determining the 
theoretical neutron star maximum mass), symmetry energies (important in determining 
the typical stellar radius and in the relative n, p, e, and v e abundances) and specific heats 
(important in determining the local temperatures). These characteristics play important 
roles in determining the matter's composition, in particular the possible presence of strange 
components (such as hyperons, a kaon condensate, or quark matter). These characteristics 
also significantly affect calculated neutrino opacities and diffusion time scales. 

Immediately following the core bounce of a massive star, the PNS shrinks from more 
than 150 km in radius to less than 20 km as neutrino emission from the low-density outer 
few tenths of a solar mass of the star robs this region of pressure support. During this same 
period, from 0.1 s to 0.5 s, substantial accretion through the supernova's shock also occurs, 
which adds mass and contributes significantly to the total neutrino emission. Nevertheless, 
both this stage and the much longer Kelvin-Helmholtz evolution which follows nearly fulfill 
the condition of hydrostatical equilibrium (Burrows & Lattimer 1986; henceforth BL in the 
text). The Kelvin-Helmholtz phase, which can last for several tens of seconds, consists of 
two major evolutionary stages. First, there is a deleptonization stage in which the excess 
trapped electron neutrinos diffuse from the central regions outward through the star. The 
diffusing neutrinos preferentially heat the stellar core while decreasing the net lepton and 
proton fractions. This stage is followed by an overall cooling stage in which the entropy in 
the star steadily decreases. 

Since the work of BL, several authors (Burrows 1988; 1990, Wilson & Mayle 1989, 
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Suzuki 1989, and Keil & Janka 1995) have performed numerical simulations of the cooling 
of a PNS. Our objective here is to improve upon these calculations by incorporating recent 
developments in the EOS of dense matter and neutrino interactions, and to perform a 
systematic study of the effects of varying the neutron star mass and the characteristics of 
the EOS. For this study, we will use the diffusion approximation for neutrino transport both 
because it provides a comparison to earlier works and because it allows us to easily assess 
how global characteristics of the neutrino emission such as average energies, integrated 
fluxes, and time scales change in response to various input parameters. This approach is 
sufficient to allow us to establish the connections between the microphysical ingredients 
and the duration of the deleptonization and cooling time scales, and to estimate the effects 
on neutrino signals. However, detailed simulations using a multigroup approach and more 
exact transport near the stellar surface will be necessary to calculate the emergent neutrino 
spectra precisely enough to make detailed predictions for current and future neutrino 
detectors. Such calculations are in progress. 

Compared to the earlier works quoted above, several improvements are carried out in 
the neutrino transport scheme. These include 

(1) the development of the appropriate diffusion equations and a derivation of explicit 
expressions for the fluxes in terms of energy, temperature-, density-, and composition- 
dependent diffusion coefficients, which can be obtained from the reaction rates of neutrinos 
in matter. We present the calculations of these coefficients using the results of Reddy, 
Prakash, & Lattimer (1998) (henceforth RPL). 

(2) the development of a transport scheme that uses the appropriate energy-averaged 
diffusion coefficients. 

The most important improvements in the EOS are that 

(1) the effects of finite entropy are calculated without approximation, and 

(2) the influence of the possible presence of strangeness in the form of hyperons is explored 
consistently. 

We employ EOSs with more realistic compressibilities and high-density symmetry energies 
and nucleon effective masses as compared to BL. Keil & Janka (1995) (henceforth KJ) 
considered the consequences of having hyperons appear in the star, but with an approximate 
treatment of the effects of finite entropy and no modification to the opacities that reflected 
the presence of hyperons. 

Several improvements over existing works are also made both in the calculation and in 
the implementation of neutrino opacities. The opacities are governed by charged current 
processes (absorption on baryons) and neutral current processes (scattering off baryons and 
electrons) in dense matter. The most significant of these improvements are: 
(1) an exact treatment of the kinematics of the neutrino-baryon and neutrino-electron 
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processes including the full effects of Pauli-blocking for arbitrary matter degeneracy, and 
(2) calculations of the effects of strong interactions arising from mass and energy shifts in a 
multicomponent system (including the presence of hyperons). 

We have found that these improvements each tend to increase neutrino mean free paths 
(RPL), and, consequently, produce shorter emission time scales relative to the situation in 
which they are ignored. 

In both BL and KJ, opacities were interpolated between expressions valid for a mixture 
of free Fermi gases in the degenerate and nondegenerate regimes, respectively. However, 
RPL demonstrated that interpolation can result in significant deviations from the exact 
phase-space integrals. RPL also showed that the exact phase-space integrals are relatively 
easy to compute, which obviates the need for interpolations. 

The calculations of Suzuki (1989) employed free gas opacities from Bruenn (1985) which 
are valid when the matter is essentially non-degenerate and non-interacting. However, for 
the conditions of interest, neither of these conditions are generally satisfied. As discussed in 
RPL, this generally results in an underestimate of the neutrino mean free paths. 

Furthermore, to include the effects of strong interactions, previous workers employed 
a constant, i.e., density-, temperature-, and composition-independent, reduction factor 
multiplying the opacity. The effects of degeneracy and strong interactions on the opacities 
were not calculated consistently with the assumed underlying EOS. 

In this paper, we focus on simulations of constant baryon mass and spherically 
symmetric protoneutron stars. We explicitly ignore effects due to accretion, rotation, 
magnetic fields and convection, which, in conjunction with multidimensional neutrino 
transport, may play an important role in the early evolution of the cores of protoneutron 
stars. In many current supernova models (Burrows 1988, Burrows & Goshy 1993), most 
accretion is completed within 0.5 s. In the constant baryon mass models we present here, it 
is assumed that the initial protoneutron star has a high-entropy mantle, so that the outer 
part of the star is initially at large radius. We find that it takes a few tenths of a second for 
the neutrinos in the high-entropy mantle to leak from the star and for the mantle to collapse 
to small radii. In effect, this largely simulates the effects of accretion on the protoneutron 
star's evolution. We find that a large fraction of the star becomes connectively unstable 
after a few seconds. The omission of convection will likely lead to an underestimate of 
neutrino luminosities and average energies at these times, and should speed up protoneutron 
star cooling compared to our models. The inclusion of better treatments of accretion and 
the inclusion of convection will be studied in a later paper. 

In §2, we develop the neutrino transport equations appropriate for the PNS problem 
by using the method of moments to solve the Boltzmann transport equations. These are 
simplified by employing the diffusion approximation, which results in expressions for the 
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number and energy fluxes in terms of energy, temperature-, density-, and composition- 
dependent diffusion coefficients. In §3, we describe the EOS employed in this paper, which 
is a finite-temperature, field-theoretical model in which the interactions among baryons 
are mediated by the exchange of a, u, and p mesons. It permits us to easily vary the 
matters' compressibility, symmetry energy, and baryonic specific heats, as well as allowing 
the systematic inclusion of hyperons, a kaon condensate, and/or quark matter. In §4, the 
opacities of matter are calculated using this EOS and the diffusion constants are determined. 
Our results for PNS evolutions are detailed in §5, in which we first describe a baseline 
simulation of a 1.6 M Q baryon mass, nucleonic EOS, PNS. We proceed to examine the 
consequences of varying the entropy and lepton content of the initial model, the inclusion 
of hyperons, and variations in the initial mass and the EOS. This section also contains a 
discussion of the expected luminosities and identification of connectively unstable zones. 
Our conclusions and outlook are contained in §6. 



2. EVOLUTION EQUATIONS 

2.1. Neutrino Transport Equations 

The equations that govern the transport of energy and lepton number are obtained 
from the Boltzmann equation (BE) for massless particles given by (Lindquist 1966) 



\dxP w dp a 



coll 



where / is the invariant neutrino distribution function, p a are the components of the 
neutrino 4-momentum with respect to the coordinate basis, are the Christoffel symbols 
for the metric. (For the most part, we will use units setting U — 1, c = 1). The indices 
a, f3, and 7 take the values t,r,8, and $. The collision term on the right hand side is due 
to neutrino-matter interactions, and is most easily evaluated in a frame comoving with the 
matter. This entails rewriting the BE in terms of the 4-momentum components of the 
neutrino with respect to the matter p a instead of the components in the coordinate basis. 
The comoving basis {e a } satisfies e a ■ e& = i] a b, with e = u representing the 4-velocity of 
matter. Using the metric 

rfs 2 = _ e 24> dt 2 + e 2A dr 2 + r 2 dQ 2 + r 2 ^2 g ^2 ^ ^ 

and taking u = (7e~*, 71* e~ A , 0, 0), where v is the fluid velocity and the Lorentz factor 
7 = (1 — f 2 )" 1 ^ 2 , the BE in the comoving basis takes the form 

(3) 




-7- 



where T^ c (a, b,c = 0, 1, 2, 3) are the Ricci rotation coefficients. The non-zero Ricci 
coefficients in the spherically symmetric case are 

(d{>ye+) , d( 1V e A )\ r0 _ rl _ _ _ A ^( 7 e A ) 

10 ~~ e e qT 1 Q ) 



<9t <9r y 

f 2 _ r _ r 3 _ pO _ _ pi _ p2 _ _ pi _ p3 _ l ve A p2 _ _p3 _ Cot (a\ 
1 20 ~~ 1 22 — 1 30 ~~ 1 33 ~~ wi 22 — UL 21 — UL 33 ~~ ui 31 ~~ > 1 33 ~~ 1 32 — • l^J 

The time t refers to that measured by an observer located at infinity. In spherical symmetry, 
the components of the neutrino 4-momentum are 

p a = (uo, lo/j,, lo(1 - /i 2 ) 1/2 cos $, co(l - /i 2 ) 1/2 sin , 

where fi is the cosine of the angle the neutrino momentum makes with respect to the radial 
direction and uo is the neutrino energy measured by a comoving observer. Then, the general 
relativistic BE in the spherically symmetric case reads (Thorne 1981) 

u{4 + t>4)% + ^K + K)^-^ 2 (^o + /^n + (i-^)r 2 2 o)|^ 

- ^(i-/x 2 )(r 00 + ri 2 + / irl - / ir 2 )g=^ . (5) 

A method often used to simplify the BE is to work with its angular moments (Thorne 
1981). This is done by applying the operator 

1 

-J dfifx 1 , i = 0,1,2,- ■■ 

to equation fl5|). Defining the i th moment of the distribution function by 

1 r+i 
Mi= 2J-i rfyUyU7 ' 
the equations corresponding to the first two moments (i = and 1) are: 
/ t dM r dM t 8M 1 r dM 1 \ 



uo 2 



' - dM x _ r2 ,dM 2 2 8M 

1 00 a. . ~+~ \ L 10 1 20J a. . + 1 20" 



duo ' v 1U M/ duo duo J 

+ uo((T 1 10 -T 2 20 )(M -3M 2 )-2(T 1 00 + T 1 22 )M 1 )=Q (6) 



( t dM l r dM 1 t dM 2 r dM 2 \ 

J x dM 2 x 2 dM 3 2 SMi 

~~ W 1 00^— + I 1 10 - 1 20J^— + 1 20 



<9u; ' v iU duo ' M duo ) 

+ w (2(r} - r 20 )(Mi - 2M 3 ) + (r 00 + ry (m - 3M 2 )) = q x , (7) 
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where 

(if 
» \Tr 



1 r+1 f df\ 1 f +1 



coll 



Multiplying the O th -moment equation (|6]) successively by (u/2n 2 ) and (o; 2 /27r 2 ), and 
integrating over the neutrino energy, one obtains the energy-averaged neutrino number and 
energy transport equations: 



t dN u r dN u t dF v r dF v 
eh^r- + el— he— - + e 



dt " dr ^ dt " 1 dr 

+ (T\ + 2Tl )N u + (Tl -2Tl 2 )F u = S N , (9) 



en^- + eA-^ \-el— — + e 



dt " 3r 1 eft " 1 3r 

+ (r^ + 3i^)j v + (2rJ -2ri,)J3 I , + (r} -it D )Pv = ^. (10) 



Above, 

f 00 du %r „ f°° du %r , „ f°° du %r n „ r°° dtu „ 9 

(12) 

N u , F u , and are the number density, number flux and number source term, respectively, 
while J u , H u , P u , and S B are the neutrino energy density, energy flux, pressure, and the 
energy source term, respectively. 

For spherically symmetric situations, equations @ and ( |10D are the transport equations 
which contain full general relativistic effects including fluid velocity terms. These equations 
must be solved together with the hydrodynamical equations to study dynamical problems 
such as core collapse supernova. For the PNS problem, however, fluid velocities are small 
enough so that hydrostatic equilibrium is nearly fulfilled and the evolution is quasi-static. 

Substituting T% c in equations (|^) and (JTO) , and utilizing the continuity equation with 
the assumption of a quasi-static evolution (that is, neglecting terms of order v), we arrive at 
the well-known neutrino transport equations (BL) in a stationary metric in the form used 
so far in protoneutron star evolution: 

d(N u /n B ) d(eHnr 2 F u ) ^S N 

+ * = e — ( 13 ) 

dt da n B 

d(J v /n B ) d(l/n B ) _, d(e 2 Hnr 2 H u ) _ ,Se n 
dt dt da n B 

where n B is the baryon number density and a is the enclosed baryon number inside a sphere 
of radius r. 
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2.2. The Equilibrium Diffusion Approximation 

At high density and for temperatures above several MeV, the source terms Qo and Q\ 
in equations (^) and ([7|) are sufficiently strong to ensure that neutrinos are in thermal and 
chemical equilibrium with the ambient matter. Thus, the neutrino distribution function in 
these regions is both nearly Fermi-Dirac and isotropic. We can approximate the distribution 
function as an expansion in terms of Legendre polynomials to 0(/i), which is known as the 
diffusion approximation. Explicitly, 

f(u,fj.) = fo(u)+ t if 1 (u), f = [l + e(^)]-\ (15) 

where f is the Fermi-Dirac distribution function at equilibrium (T = T mat , \i v = /i^ 9 ), with 
lu and fi u being the neutrino energy and chemical potential, respectively. 

The main goal is to obtain a relation for fi in terms of /q. In this approximation, the 
moments Mi of / which appear in equations (^j) and ((?]) are 

M = / , M 1 = \h, M 2 = ~/ 0) and M 3 = \f ± . (16) 

6 6 5 

We substitute these into equation (|7|) and again neglect all fluid velocity dependent terms 
(i.e., eQ,e\, and T^q). Furthermore, since the transport is driven by slowly varying spatial 
gradients in temperature and chemical potential, the time derivative of f\ has also been 
neglected. One then finds 

e _ A fdfo _ UJ ^^fo\ = 3 Qi / 17 x 
\dr dr du ) u 

The interaction between neutrinos and matter is dominated by scattering (from both 
baryons and leptons), absorption (on baryons) and emission (due to electron capture) 
processes; thus the collision term Qi is given by 

(IL="(^- fl -£ +wi - fl -£)- (i8) 

where j a is the emissivity, A a is the absorptivity, and j s and A s represent the scattering 
contributions to the source term. These are 



roc rl rZ7Y 

/ du'u' 2 / du' I d$ f(u', ft') B^(u,u', cos 6), (19) 
J o J-x Jo 

T- = 7^3/ dJj2 \ d^[l-f(u , ,fj!))Ii^ t (u,u/,coae), (20) 

A s ( Z7T ) -'0 J—l JO 

where 6 is the angle between the ingoing and outgoing neutrino. From the reciprocity 
relations between a reaction rate and the rate of its inverse process (the principle of detailed 
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balance), the emissivities and absorptivities (R™ and _R° U * in the scattering kernels) are 
related according to 



A a (o;) 



a (u) and R™ = e^ u '- u3) R™ t 



(21) 



The different elementary processes which contribute to equation fllq) are described in §4. 

To arrive at the final form for the diffusion equation and the appropriate definitions for 
the diffusion constants we expand the scattering kernel in terms of its Legendre moments: 



R° ut = J dcos6Pi{cos0)R™ t {u,u',cosl 



Then, one finds 



Qo = u [j a (l - /oM) 



Qt 



-u- 



AH 



3a + ^ I" K l I • 



(22) 



(23) 



where the scattering contributions to Qo and Qi are given by 
1 



(2tt) 2 -A) 



du/f/ 1 {[/oCf/)^-") - / (w) + (1 - e^-^/oCaO/oCw) 



ryout 

K 



+ (i _ ^k-^)|M4AH^ 



out 



(24) 



(2tt) 2 

AH) 

AH 



^' W ' 2 {[l-/ ( W ')(1 



— e 



/3(o/-u/h 



ryout 

K 



*'- u) +/o(w)(l-e^ ) ) fl? 



out 



In the case that the scattering is isoenergetic, k\ may be further simplified: 

,2 ,1 

y d cos (1- cos 0)^(^,0;, cos 0), 



L) 



(27T) 



(25) 



(26) 



which is the well-known transport opacity used in many applications. If we substitute the 
above expression for Qi into equation (|T7D, we obtain the relation between f and fi that 
we are looking for: 



A 



-D(oj) 



c -A d fo ue -A^_^fo 

dr dr du 



The explicit form of the diffusion coefficient D appearing above is given by 



D(u) 



3 + K +Kl 



-1 



(27) 



(28) 
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Because in most situations of interest to us neutrino scattering is dominated by 
scattering off baryons, the scattering processes are nearly elastic and the implicit 
dependence on fi in the diffusion coefficient may be safely ignored. However, should 
electron scattering or the inelastic contribution to baryon scattering become important, one 
is forced to solve a more complex integro-differential equation than the familiar Fick's law. 

To obtain explicitly the equations that determine the lepton number and energy fluxes 
we note that 



dfo 
dr 



T ^ + ^dT\dfo 
dr T dr J doj 



(29) 



where r\ v 
equation 



fi u /T is the neutrino degeneracy parameter. Substituting the above relation in 
7|), we obtain 



fi = -D{u)e 



T dr] uj d(Te 4> ) 
dr Te^ dr 



dfo 
duj 



(30) 



In line with the definition of the fluxes given before in equations ( p~T| ) and (|T"2|) , the 

energy-integrated lepton and energy fluxes are 



e~ A e-*T 2 



Qn 2 



e~ A e-*T 3 



6n 2 



dr dr 
dr dr 



(31) 



The coefficients D 2 , D 3 , and D4 are related to the energy- dependent diffusion coefficient 
D(u) through 

D n = I dx x n D(cu)f (Lu)(l - f (u)) , (32) 



where x = u/T. Note that these diffusion coefficients arise naturally from the transport 
equations. They depend only on the microphysics of the neutrino-matter interactions (see 
§4 for details). The fluxes appearing in the above equations are for one particle species. To 
include all six neutrino types, we redefine the diffusion coefficients in equation (|3T|): 



D 2 = D l n 



D 



2 5 



D 



3 > 



D 3 = D? + D^ + ADt . 



(33) 



2.3. Numerical Notes 

We generally assume in this paper that the total baryon number in each simulation 
is conserved, that is, we do not include the effects of accretion. These effects will be 
incorporated in future work. The enclosed baryon number a is then a convenient Lagrangian 
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variable. The equations to be solved split naturally into a transport part, which has a 
strong time dependence, and a structure part, in which evolution is much slower. We choose 
to solve these two parts separately in an iterative fashion which is equivalent to solving 
them together. 

Explicitly, the structure equations are 

dr 1 dm p . ^. 

da 47rr 2 riBe A ' da riBe A 



■(m + A7rr 3 P) , — = -(p + P) (m + 47rr 3 P). (35) 



da ATcr 4 nB^ ' da Airr^riB 

The quantities m, p, and P include contributions from the leptons. To obtain the equations 
employed in the transport, we combine equation and the corresponding equation for 
the electron fraction 

^ = (36) 
dt n B 

to obtain 

«£ +e -*^W) =0 . (37) 
at da 



Similarly, we combine equation (O) with the matter energy equation 



dU + p d(l/n B ) = _ e ,S_B 
at dt ub 

where U is the specific internal energy and use the first law of thermodynamics to obtain 

In regions of small optical depth, diffusion codes can predict fluxes that exceed the 
black-body limit. To prevent this, we correct the fluxes by a flux limiter 3A(£), where 
= F u /N u and = H u /J u for the number and energy fluxes, respectively. For the 
functional form of A(£), we tested both the Levermore & Pomraning (1981) and the Bowers 
& Wilson (1982) flux-limiters. Only the outermost few mass shells are affected by this 
choice, and we have not found it to make an appreciable difference in the overall evolution or 
in the mean neutrino energies or luminosities. For the PNS cooling problem, as long as most 
of the transport is nearly diffusive, small errors in the transport near the surface cannot 
affect the global evolution. This, however, is not the case in supernova simulations, in 
which the flux-limiter affects more of the computational grid and an appropriate treatment 
becomes a crucial aspect of transport. 
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The structure equations ( p4|) and (|35| ) satisfy the boundary conditions 



r( a = 0) = 0; m(a = 0) = 0, 




P(a = a,) = P t 



S ) 



(40) 



where a s is the total number of baryons in the star, which remains fixed during the 
evolution, and P s is the surface pressure, which is chosen to be very small. We have checked 



the center (this is equivalent to imposing null gradients of the temperature and chemical 
potential). At the surface, we may write 



where a n and a e are constants chosen to each be in the range 0.4 — 0.6. The results for 
overall neutrino number and energy fluxes are insensitive to the exact choice of these 
constants, due to the fact that the fluxes saturate near the surface and the values of Y v in 
the surface region adjusts itself to compensate for small changes in a n and a e . 

An initial model is specified by choosing entropy s and lepton fraction Yl profiles in 
a and relaxing the star to hydrostatic equilibrium, which determines values for the radius 
r, gravitational mass m, gravitational potential 0, and pressure P. In general, this initial 
model will not be in radiative equilibrium; however, we found that radiative equilibrium is 
established within about 5 ms, a time too short to be important for the protoneutron star 
evolution. 

The evolution is begun by taking an initial transport time step, determining the 
quantities Y[(a) and s'(a) as first estimates for the lepton fraction and entropy at the 
end of the time step. In this step, the radius, gravitational mass and pressure profiles are 
kept fixed. Now using these values of lepton number Y[(a) and entropy s'(a), we solve 
the structural equations relaxing the star into hydrostatic equilibrium, which determines 
r'(a), m'(a), <f) f {a), and P'(a). Using these primed values for r, m,0, and P, but the original 
values of s(a) and Yl(o), we retake the transport time step, thereby determining updated 
estimates, s"(a) and Y^(a). This is followed by another hydrostatic equilibrium relaxation 
which determines r"(a), m"(a), <f)"{a), and P"(a). This iterative procedure is found to 
converge rapidly, typically in 3-4 iterations. This method is equivalent to using a full 
implicit discretization scheme to solve the structure and transport equations simultaneously. 




£f = a n and = a, 



(41) 
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2.4. Neutrino Luminosities 

The ultimate objective of an observation of neutrinos from a Galactic supernova will 
be to deduce the properties of the PNS, including its mass, composition, and the EOS. The 
prediction of neutrino signals, and the determination of how well the mass and properties 
of the EOS can be estimated from them, is an important goal of our work. The signal 
depends both upon the neutrino fluxes and energy spectra. Its accurate determination 
requires the use of better transport than the diffusive, approach used in this paper. 
Thus, we have limited our predictions to those utilizing information only from the total 
neutrino luminostiy and the average emitted energy of all neutrinos. However, our neutrino 
transport calculations tracked each flavor of neutrino and the diffusion constants were 
defined appropriately. In fact, detailed transport calculations show that the emitted spectra 
are, to a fair approximation, Fermi-Dirac with a small effective degeneracy parameter for 
all flavors of neutrinos (Mayle, Wilson, & Schramm 1987). Furthermore, the presence of a 
small degeneracy parameter gives negligible corrections to a detector's count rate (Lattimer 
& Yahil 1989). Thus, a fair representation of the signal in a terrestrial detector can be found 
from the time dependence of the total neutrino luminosity and average neutrino energy 
together with an assumption of a Fermi-Dirac spectrum with zero chemical potential. 

The total neutrino luminosity is the time rate of change of the star's gravitational 
mass, and is therefore primarily a global property of the evolution. This luminosity, due to 
energy conservation, must also equal 

L v = e 2 H7rr 2 H u (42) 

at the edge of the star. We have checked this test of energy conservation for our code 
and find that it is valid to better than 1% at all times greater than about 5 ms, when 
the star comes into radiative equilibrium. For times greater than about 5 ms, initial 
transients become quite small and the predicted luminosities should be relatively accurate 
compared to full transport simulation. We estimate the average energy of neutrinos from 
the temperature T v of the matter at the neutrinosphere R u , defined to be the location in the 
star where the flux factor = 0.25. However, since the spectrum may not be Fermi-Dirac 
at the neutrinosphere, our transport scheme cannot give a very precise value for the average 
energy. Comparing with better transport in this region, we find that the average energy 
< E v >Ri 3T U , where T v is a mass average in the outermost zone. Because it is a globally 
determined quantity, the luminosity L v is necessarily more accurately determined than 
either R v or T v . In lieu of more complete transport results, we will place more weight on 
L v in the estimate of detector signals in §5.7. 
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3. THE EQUATION OF STATE 

The role of the EOS in determining the structure of neutron stars has been studied 
extensively (see, for example, Prakash et al. 1997 & references therein). The masses 
and radii of neutron stars depend upon the matters' compressibility, the composition of 
matter at high density, and the nuclear symmetry energy. In the PNS problem, the finite 
temperature aspects of the EOS also play an important role. During the early evolution 
the entropy in the central regions is moderately high s ~ 1 — 2, which correspond to 
temperatures in the range T = 20 — 50 MeV. In this work we employ a finite temperature 
field-theoretical model where the interactions between baryons are mediated by the 
exchange of cr,u, and p mesons. Including the leptonic contributions, the total Lagrangian 
density is given by (Serot & Walecka 1992), 

L = Li + Ln 

i i 
- \w^W^ + \mlu^ - ^B^B"" + \m 2 p b^ + \d^a - ^m 2 a a 2 - E%)(43) 

Here, the leptons, / = e~ and /i - , and B are the Dirac spinors for baryons and t is the 
isospin operator. The sums include baryons i = n,p,A, E, and S. The field strength tensors 
for the u) and p mesons are W^ u = d^u u — d u u p and B Mly = d^b u — d u b^, respectively. The 
potential U (a) represents the self-interactions of the scalar field and is taken to be of the 
form 

U(a) = h)M n {g aN af + ^c(g aN a) 4 . (44) 

Electrons and muons are included in the model as noninteracting particles, since their 
interactions give small contributions compared to those of their free Fermi gas parts. 

In the mean field approximation, the partition function (denoted by Zjj) for the 
hadronic degrees of freedom is given by 

In Z H = PV [|m> 2 + - \™> 2 - U{a)\ + 2V £ / ^ In (l + e^ E ^) , (45) 

where (3 = T~ l and V is the volume (we take ks = 1)- The contribution of antibaryons is 
not significant for the thermodynamics of interest here, but is straightforwardly included in 



equation (f45|) . Here, the effective baryon masses M* = Mi — g a i<j and E* = yk 2 + M* 2 . 
The chemical potentials are given by 



Pi — Vi + g W iU)Q + gpitsibo , 



(46) 
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where t 3i is the third component of isospin for the baryon. Note that particles with t 3i = 0, 
such as the A and S° do not couple to the p. The effective chemical potential V{ sets the 
scale of the temperature dependence of the thermodynamical functions. 

Using Zh-i the thermodynamic quantities can be obtained in the standard way. The 
pressure Pjj = TV' 1 liaZ H , the number density for species B, and the energy density e# 
are given by 



71: 



d3k ( e m-»i) + iY 1 



en = \mW + U(a) + §m> 2 + frffi + 2 £ / (e^~^ + l)^ . (47) 

The entropy density is then given by sh = (3{£h + Ph — J2if^i n i)- The meson fields are 
obtained by extremization of the partition function, which yields the equations 

m ^o = J2 9^ ; 7n 2 p b = J2 QpihBrii , 



2__ du i°) , \- . . r d3k MlLm-»i) ■ ir 1 



The total partition function Z tota [ = ZhZ^, where Z^ is the standard noninteracting 
partition function of the leptons. The additional conditions needed to obtain a solution are 
provided by the charge neutrality requirement, and, when neutrinos are not trapped, the 
set of equilibrium chemical potential relations required by the general condition 

Ik = kUn - qiHi , (49) 

where hi is the baryon number of particle i and is its charge. For example, when i = e~, 
this implies the equalities 

A*a = = /W S o = fi n , u E - = fi B - = u„ + fi e , and fi p = « s + = fi n - fi e . (50) 

In the case that the neutrinos are trapped, equation (§^) is replaced by 

Hi = bifin - qi(fj,i - ii Vi ) . (51) 

The new equalities are then obtained by the replacement \i e — > « e — fi Ue in equation 
(|50"D. The introduction of additional variables, the neutrino chemical potentials, requires 
additional constraints, which we supply by fixing the lepton fractions, Y Li , appropriate for 
conditions prevailing in the evolution of the protoneutron star. The contribution to pressure 
from neutrinos of a given species is given by their free Fermi gas expressions. 



In the nucleon sector, the constants g a ^, g^N, 9 P N>b, and c are determined by 
reproducing the nuclear matter equilibrium density no, and the binding energy per nucleon, 
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the symmetry energy, the compression modulus K 0l and the nucleon Dirac effective mass 
M* at uq. The parameters of the models are listed in Table 1 and for the most part are 
taken from Glendenning & Moszkowski (1991). Model GM1 is a stiff EOS (K = 300 MeV) 
while the model GM3 is softer with K = 240 MeV. To study the possible role played by 
the nuclear symmetry energy we also investigate the model GM4, which is similar to GM3 
but with a smaller symmetry energy. 

In addition to models containing only nucleonic degrees of freedom (GMlnp & GM3np) 
we investigate models that allow for the presence of hyperons (GMlnpH & GM3npH). 
The hyperon coupling constants may be determined by reproducing the binding energy of 
the A hyperon in nuclear matter (Glendenning & Moszkowski 1991). Parameterizing the 
hyperon-meson couplings in terms of nucleon-meson couplings through 

x aH = 9aul 9aNi %uH = QuH 1 9uNi x pH = QpUj ' 9pN , (52) 

the A binding energy at nuclear density is given by 

(B/A) A = -28 = x uA g uN u - x aA g aN a , (53) 

in units of MeV. Thus, a particular choice of x a \ determines x^a uniquely. To keep the 
number of parameters small, the coupling constant ratios for all the different hyperons are 
assumed to be the same. In model GM1, the hyperon couplings are given by 

x a = XaA = XoY, = x aZ = 0.6 , and X u = X uA = XuY, = x^-e = 0.653 . (54) 

The p-coupling is of less consequence and is taken to be of similar order, i.e. x p = x a . In 
GM3, x (j = 0.6 , x u = 0.659 and x p = 0.6. 

In Figure [I] the particle fractions for models GM3np and GM3npH are shown with 
and without trapped neutrinos. The left panels show the particle fractions in hot and 
neutrino-free matter (mimics the ambient conditions subsequent to deleptonization) and 
the right panels show the particle fractions in neutrino trapped matter (ambient conditions 
in the core at birth). Considering matter with nucleons-only first, note that the major effect 
of trapping is to make the electron (and hence the proton) concentration high compared to 
the case in which matter is neutrino free. Compared to the free gas case the proton fraction 
is appreciably larger in matter with strong interactions due to the larger symmetry energies. 
At fixed Y L , this acts to decrease the electron neutrino fraction at high density (see top right 
panel). During the deleptonization phase, the EOS affects the evolution mainly through 
the density dependent nuclear symmetry energy as it directly influences A„ e (see §4) and 
the magnitude of the spatial gradient dfj, u /dr (see §2). Subsequent to deleptonization, the 
proton fractions decreases, but is still high enough to enable the charged current reactions 
to occur for most relevant densities and temperatures of interest (see §4). 
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Although the threshold densities for the appearance of hyperons are sensitive to 
the hyperon coupling constants (Knorren, Prakash, & Ellis 1997), the higher entropies 
encountered during the early evolution generally results in hyperon fractions which are 
larger compared to the case of zero temperature matter (bottom left panel). In lepton-rich 
matter (bottom right panel), the appearance of hyperons is delayed to higher densities 
due to the relative dispositions of the various chemical potentials. With the appearance 
of the hyperons, the neutrino population begins to increase with density, in contrast to 
the monotonic decrease exhibited in the hyperon free case. Since the negatively charged 
hyperons effectively replace the electrons at high density, the final lepton content of stars 
with hyperons is smaller relative to those in nucleons-only matter. Correspondingly, we 
should expect to see a larger lepton number radiated from stars with hyperons (see §5). 

In Figure the variation of pressure with baryon density is shown. The upper panels 
show the results for s = 2 and the lower panels for s = 1. Thermal effects contribute only 
modestly to the total pressure. The pressure is dominated from effects of interactions and 
degeneracy at high density. The left panels in Figure |2| show results for Y Ue = 0, and the 
right panels for a fixed lepton fraction of Yl = 0.4. While the differences between lepton 
rich and lepton poor matter for the np models are negligible, significant differences arise 
in models with hyperons. When the electron neutrino fraction is large, their presence 
inhibits the appearance of hyperons (see also Figure [I]; note however, that a large hyperonic 
component at high density in the neutrino free case significantly softens the EOS (Prakash 
et al. 1997)). This softening with decreasing lepton number will lead to compression during 
the deleptonization epoch. Hence, there exists a window of initial masses for which the star 
becomes unstable to gravitational collapse during deleptonization. The maximum masses 
supported at low temperature and for lepton poor matter are significantly smaller than for 
hot and lepton rich matter. 

Finite temperature properties of matter at high density influence the diffusion of 
neutrinos especially through the specific heat. Neutrino mean free paths are strongly 
temperature dependent (see §4) and the ambient matter temperature controls the diffusion 
of neutrinos to a large extent. In Figure |3|, we show the matter temperature at fixed entropy 
per baryon as a function of the baryon density to contrast the behavior observed in models 
with and without hyperons. The trends observed in the np models are similar while the 
differences between np and npH models are very significant. The appearance of additional 
fermionic degrees of freedom in the form of hyperons (we may expect similar behavior in 
models with quarks (Prakash, Cooke, & Lattimer 1995)) increases the low temperature 
specific heat. This generally favors lower temperatures during the evolution for the npH 
models, all else being equal. This is because, in charge neutral and beta-equilibrated 
matter, bulk of the entropy resides in the baryons (Prakash et al. 1997). Under degenerate 
situations {T jEp i » 1), this result may be understood from the relation that connects the 



temperature and baryonic entropy (or specific heat) through the concentrations (Yj = rii/ri): 



(37r 2 n B ) 2 / 3 



for M* » k Fi 
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(3yr 2 n B ) 2 / 3 




(55) 



for M* « k Fi 



where k Fi = (37r 2 7iBYj) 2//3 are the Fermi momenta. For the temperatures of interest 
here, and particularly with increasing density, the above relation provides an accurate 
representation of the exact results for entropies per baryon even up to s = 2. The behavior 
with density of both the concentrations and effective masses controls the temperatures for a 
fixed s. In the absence of any variation with M*, a system with more components at a given 
baryon density has a smaller temperature than a system with fewer components (recall that 
J2i Yi — 1). Effective masses that rapidly drop with density oppose this behavior and lead to 
larger temperatures. The fact that the temperatures are smaller in the presence of hyperons 
than those in nucleons-only matter attests to the predominant effect of concentrations. 
Note that the differences between the temperatures for the cases with and without hyperons 
are larger than those obtained by using different models for the individual cases. 

However, as we have seen from Figure |2|, the appearance of hyperons also leads to 
higher central densities due to softening; this acts to compensate for the higher specific heat 
in npH models. In addition, the larger neutrino chemical potentials (see Figure ^) in npH 
models generate larger entropies during the deleptonization epoch; this also acts to increase 
the central temperature in these models. Results discussed in §5 clearly show that despite 
the significant differences seen in Figure |3] the maximum central temperatures reached in 
the both np and npH models are nearly equal due to the above mentioned feedbacks. 

The lepton chemical potentials influence the deleptonization epoch (see §5) and are 
shown in Figure |j. For np models a lower nuclear symmetry energy favors a larger v e 
fraction and has little effect on the e~ fraction at Y/ = 0.4. Models with hyperons lead 
to significantly larger \x Ve and lower /i e , both of which influence the diffusion of electron 
neutrinos. The electron chemical potentials in neutrino free matter are reduced to a 
greater extent by changes in composition and symmetry energy as there are no neutrinos to 
compensate for changes in (x = fi n — /i p . The impact these modifications will have on the 
evolution are discussed in §5. 

We turn now to some aspects of the global energetics. The binding energy, which is 
the difference between baryonic and gravitational masses, is an important observational 
parameter because nearly 99% of it appears as radiated neutrino energy. The rate of 
change of binding is essentially the total neutrino luminosity. In Figure ^|, we display the 
binding energy as a function the baryonic mass for the models considered here. As found 
by Lattimer & Yahil (1989) and Prakash et al. (1997), there exists a universal empirical 
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relationship between the binding energy and the mass: 

( M R \ 2 

B.E. = (M B - M G )c 2 ~ (0.065 ± 0.01) — - M . (56) 

This universality is not altered by the presence of significant softening due to the presence 
of hyperons, kaons or quarks. Since nearly all of the binding energy is released in the 
form of neutrinos, an accurate measurement of the total neutrino energy will lead to a good 
estimate of the remnant mass. However, as the results of Figure [5] show, it will not be 
possible to distinguish the various EOSs from the total binding energy alone. 

Table 2 shows the maximum gravitational masses for the models studied. The three 
columns roughly depict the three important stages of the evolution of a protoneutron star. 
At birth the PNS has s = 1 and Yl = 0.4, after deleptonization the entropy has reached its 
maximum (s ~ 2), and finally cools down to its final cold catalyzed state. 



4. NEUTRINO-MATTER INTERACTION RATES 

One of the important microphysical inputs in PNS simulations is the neutrino opacity 
at supra-nuclear density (Bruenn 1985; Mezacappa Sz Bruenn 1993; Burrows & Lattimer 
1986; Wilson & Mayle 1989; Suzuki & Sato 1992; Keil 1994; Keil & Janka 1995). However, 
calculations of neutrino opacities in dense matter have received relatively little attention 
compared to other physical inputs such as the EOS. Although it was realized over a decade 
ago that the effects due to degeneracy and strong interactions significantly alter the neutrino 
mean free paths, it is only recently that detailed calculations have become available (Reddy 
& Prakash 1997; RPL 1998; Reddy et al. 1997, Prakash & Reddy 1997, Burrows & Sawyer, 
1998a,b). The weak interaction rates in hot and dense matter are modified due to many 
in-medium effects. The most important of these are: 

(1) Composition: The neutrino mean free paths depend sensitively on the composition 
which is sensitive to the nature of strong interactions. First, the different degeneracies of 
the different Fermions determines the single-pair response due to Pauli blocking. Second, 
neutrinos couple differently to different baryonic species; consequently, the net rates will 
depend on the individual concentrations. 

(2) In-medium dispersion relations: At high density, the single-particle spectra are 
significantly modified from their noninteracting forms due to effects of strong interactions. 
Interacting matter features smaller effective baryon masses and energy shifts relative to 
non-interacting matter. 

(3) Correlations: Repulsive particle-hole interactions and Coulomb interactions generally 
result in a screened dielectric response and also lead to collective excitations in matter. 
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These effects may be calculated using the Random Phase Approximation (RPA), in which 
ring diagrams are summed to all orders. Model calculations (Sawyer 1975, f989; Iwamoto 
& Pethick 1982; Horowitz & Wehrberger 1991a,b; 1992, Raffelt & Seckel 1995, Sigl 1995, 
Reddy et al. 1997, Prakash & Reddy 1997, Burrows & Sawyer 1998a, b) indicate that at 
high density the neutrino cross sections are suppressed relative to the case in which these 
effects are ignored. In addition, these correlations enhance the average energy transfer 
in neutrino-nucleon collisions. Improvements in determining the many-body dynamic 
form factor and assessing the role of particle-particle interactions in dense matter at 
finite temperature are necessary before the full effects of many-body correlations may be 
ascertained. 

(4) Axial charge renormalization: In dense matter, the axial charge of the baryons is 
renormalized (Wilkinson 1973, Rho 1974, Brown & Rho 1991), which alters the neutrino- 
baryon couplings from their vacuum values. Since the axial contribution to the scattering 
and absorption reactions is typically three times larger than the vector contributions, small 
changes in the axial vector coupling constants significantly affect the cross sections. 

(5) Multi-Pair excitations: Neutrinos can also excite many-particle states in an interacting 
system, inverse bremsstrahlung being an example of a two-particle excitation (Hannestad & 
Raffelt 1998). These excitations provide an efficient means of transferring energy between 
the neutrinos and baryons which are potentially significant in low-density matter. However, 
multigroup neutrino transport will be needed to fully include this affect. In addition, 
such calculations require source terms for neutrino processes such as bremsstrahlung and 
neutrino pair production. The latter process has been accurately treated by Pons, Miralles, 
& Ibahez (1998). 

The relative importance of the various effects described above on neutrino transport is 
only beginning to be studied systematically. As a first step, we will focus on effects due to 
modifications (1) and (2) above, since RPL have shown how to calculate them consistently 
with the underlying EOS. Work on the other modifications listed above is currently in 
progress and will be incorporated in a future publication. In what follows, we recapitulate 
some of the essential results for the evaluation of the neutrino mean free paths and diffusion 
constants. 

The scattering and absorption reactions that contribute to the neutrino opacity are 

v e + B -> e~ + B', u e + B^e + + B', (57) 
v x + B -> v x + B', v x + e~ -> v x + e~ , (58) 

where the scattering reactions are common to all neutrino species and the dominant source 
of opacity for the electron neutrinos is due to the charged reaction. The total cross section 
per unit volume for the scattering reaction v + B — > v' + B' and absorption reaction 
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v + B —»■ e + B' are given by 

— ~ 2^ 7-oc d% 1 - exp(-z) io| ^ g g - [ ^ + ^ 2 + ^ • (59) 

where .E and E' are incoming and outgoing lepton energies, and g M = (qo,q) is the four 
momentum transfer to the baryons so that q 2 = 2EE'(cos6 — 1), where 9 is the angle 
between the incoming and outgoing lepton. The factor (1 — exp(— z)) arises due to the 
principle of detailed balance, and z = (go + — Hb>)/T. The response functions Ri, R2, 
and i?3 are functions of the energy and momentum transfers q$ = E — E' and q = k — k', and 
strongly depend on the properties of matter, such as the individual particle degeneracies 
and in-medium modifications to the baryon propagators. The density and temperature 
dependent baryon effective masses and chemical potentials required for their evaluation are 
obtained directly from the EOS. Explicitly, the response functions are written in terms of 
the neutrino baryon coupling constants V and A, and the polarization functions 

R x = (V 2 + A 2 ) [Imn£(g ,g) + Imn£(g ,g)] (60) 
R 2 = (V 2 + A 2 )lmU^(q ,q)-A 2 lmUi(q ,q) (61) 
R 3 = 2VAlmU^ A (q ,q). (62) 

In the case of absorption V = Cgy and A = CgA, where C is the Cabibbo factor, gy and g& 
are the vector and axial-vector couplings. For scattering reactions V = 5v/2 and A = gA/2, 
respectively. Numerical values are taken from Tables I and II of RPL. The imaginary part 
of the various polarization functions IT may be found in (Horowitz & Wehrberger 1991, 



1992; RPL 1998). Comparing with equation ( |20| ) we identify 



R out = 4G| (cosg [ARi + R 2 + BR 3 ) (63) 



I ( t r°° (ens — 1 1 

= _£ / dEe El [1 - f eq (E e )] / d cos 6 i- -i + i? 2 + 5i? 3 ] (64) 

A a 7r jo j-i 1 — e 

Note that the above formalism accounts for the effects arising due to matter degeneracy 
and in-medium modifications in the baryon dispersion relations exactly. As stated earlier, 
ignoring the effects due to correlations, the total scattering kernel R° ut is simply the 
incoherent sum over all possible scattering process. Similarly, the absorption mean free 
path A a is the sum over all possible absorption reactions. 

Under degenerate conditions even modest changes to the composition significantly alter 
the neutrino scattering and absorption mean free paths. In Figure the neutrino scattering 
and absorption cross sections are shown for models GM3np and GM3npH relevant to the 
deleptonization and cooling epochs. The top panels show the scattering cross sections 



-23- 



common to all neutrino species in neutrino free matter. The scattering cross section per unit 
volume for thermal neutrinos (E u = ttT) is shown in the left panel for various temperatures. 
To study the influence of hyperons, the ratio of the cr npH /a np is shown in the right panels. 
The presence of hyperons significantly increase the scattering cross sections, by a factor 
~ (2 — 3). Similar results for the absorption cross sections are shown in the lower panels for 
Yl = 0.4. Again we notice a significant enhancement (right panel) when hyperons appear, 
the factor here could be as large as 5. 

During the deleptonization stage, lepton number transport is sensitive to charged 
current reactions which dominate scattering reactions. At zero temperature, charged 
current reactions v + n «-> e + p depend sensitively on the proton fraction Y p (Lattimer et al. 
1991). Kinematic restrictions require Y p to be larger than 11 — 14% (direct Urea threshold). 
At early times, a finite neutrino chemical potential favors a large Y p throughout the star, 
which enables these reactions to proceed without any hindrance. Toward the late stages, 
however, Y p decreases with decreasing p> v and charged current reactions may be naively 
expected to become inoperative. The threshold density for the charged current reaction 
when fi u = and T = depends sensitively on the density dependence of the nuclear 
symmetry energy. In field-theoretical models, in which the symmetry energy is largely given 
by contributions due to p-meson exchange, the critical density is typically n B = 2 ~ 3n . 
However, finite temperatures favor larger Y p s and increase the average neutrino energy 
enabling the charged current reactions to proceed even below these densities. Figure |7| 
shows that this is the case even at relatively low temperatures (T ~ 3 — 5) MeV for a 
baryon density n# = 0.15 fm~ 3 . The sharp rise with temperature, which occurs even for 
Y v = 0, clearly indicates that this reaction dominates the v e opacity even during the late 
deleptonization era. Thus, charged current reactions cannot be simply turned off when the 
neutrino chemical potential becomes small enough as was done in prior PNS simulations 
(BL). 

The EOS and neutrino mean free paths are intimately related, which is best illustrated 
by comparing the results shown in Figure |6| with those shown in Figure [8]. Composition and 
the baryon effective masses influence both the neutrino mean free paths and the matter's 
specific heat. Hyperons decrease the neutrino mean free paths at constant temperature 
(Figure ||). This trend is reversed at constant entropy due to the significantly lower 
temperatures favored in npH matter. Similar effects are apparent when we compare np 
models with different baryon effective masses. At a constant temperature, the larger 
effective mass in model GM3np favors larger cross sections, while at constant entropy this 
trend is again reversed due to the lower temperatures favored by the larger specific heat. 

The diffusion coefficients are calculated using equation (|32|) with the cross sections 
discussed above. Earlier works (BL Sz KJ) made simplifying assumptions to define the 
diffusion coefficients including the use of interpolation to bridge the degenerate and 
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nondegenerate limits. But in partially degenerate matter, effects due to Fermi motion and 
Pauli blocking are substantial. The effects of density-dependent baryon effective masses, 
which are significant, were also ignored. In this work, we numerically integrate equation fl3~2|) 
over the neutrino energies by employing the exact cross sections. The diffusion coefficients 
D2,D%, and D4 are functions of Ub, T, and Y Ue . 

5. RESULTS AND DISCUSSION 

5.1. Baseline Model 

We have performed simulations for several dense matter models, initial conditions and 
baryonic masses. To begin, we describe in detail the general features of a reference model, 
corresponding to the dense matter EOS containing nucleons and leptons (model GM3np), 
and for a baryon mass Mb = 1.6 M . We shall initially concentrate on the temporal 
evolution of the various thermodynamic quantities, such as s, T, ns, Y e , Y v and fi u inside 
the PNS (see Figure |9|). Our overall results are qualitatively similar to those obtained in 
earlier works (BL; KJ), but some interesting differences exist nevertheless. 

The initial entropy and Yl profiles were taken from the profiles obtained in detailed core 
collapse calculations of an Mb = 1-08 M Q star (Wilson & Mayle 1989) within a hundredth 
of a second of core bounce. For larger mass PNSs, the s and Y L profiles were scaled in 
terms of the relative enclosed baryon mass. At early times Yl ~ 0.36 and s ~ 1 in the 
central regions, while the shocked outer mantle which begins at about 0.5 M (for the 1.08 
M Q star) is at high entropy (s > 8) and contains smaller lepton fractions. The outermost 
regions, which contains subnuclear density matter, and tends to nearly symmetric matter at 
low density, hence, it contains a large positive gradient in electron concentration. By scaling 
the results for 1.08 M for larger masses, we attempted to simulate the profiles expected for 
larger collapsing stars. The masses of both the inner core with roughly constant s and Yl 
and the shocked, relatively lepton-poor outer mantle should increase with total mass. 

The evolution of the 1.6 M star proceeds as follows. The steep neutrino 
chemical potential gradients drive a radially outward flow of electron neutrinos, causing 
deleptonization. The most rapid deleptonization occurs in the low-density outer mantle, 
which collapses during the first few tenths of a second due to substantial pressure loss from 
neutrino emission. Note from Figure |^ that compression of the outer half of the star leads 
to a significant temperature increase. The deleptonization of the central regions takes much 
longer, tens of seconds in this example. The central regions, as they deleptonize, are also 
heated due to neutrino transport (Joule heating). Eventually, after about 15 seconds, the 
central heating is overwhelmed by neutrino cooling and the central temperature begins a 
steady decrease. This star took about 50 seconds to cool to below 5 MeV temperatures in 
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the central regions. At these low temperatures, the neutrino mean free path has increased to 
values approaching that of the stellar radius and the star became transparent to neutrinos. 
Beyond this time, diffusion is no longer a valid approximation to the transport equations; 
however, the neutrino luminosities are very low (see below) and perhaps even undetectable. 

Qualitatively, the simulations of BL are similar: the central entropy rises from s = 1 
to the peak value of about s = 2 in about 15 seconds, and neutrino transparency occurs 
after about 50 s. However, one significant difference is that at the beginning of the central 
cooling epoch, when the central temperature reaches a maximum, there are still substantial 
numbers of trapped neutrinos present (fi u » T). Consequently, in our simulations the 
deleptonization continues through most of the cooling epoch as well. We attribute this 
difference to differences in the treatment of the nuclear absorption processes. In the 
calculations of BL, the absorption cross section does not contain any explicit temperature 
dependence and, further, is assumed to vanish when Y e fell below a critical value of about 
1/9. The value 1/9 corresponds to the critical proton fraction below which the direct Urea 
process is prohibited at zero temperature (Lattimer et al. 1991), in the absence of muons 
(both BL and this paper neglect muons). However, as shown in Figure [7], the absorption 
cross section does not behave as a step function as long as the temperature is above a few 
MeV. In addition, BL used a symmetry energy with almost no density dependence, with the 
result that fi n — fi p , and consequently the equilibrium value of the proton fraction, decreased 
to the critical value relatively early in the deleptonization. Hence, in our calculation, 
absorption contributes to the opacities for a much longer time than in BL, so that the 
decrease of \i v in the central regions is suppressed. 

Despite the differences in the evolution of /j, u , the long-timescale cooling is dominated 



by mu- and tau-neutrinos. The upper panels of Figure |10| shows neutrino mean free paths 
for an average energy, < E u >= jj, v + nT, for four times during the evolution. The mu- and 
tau-neutrino mean free paths are significantly larger than those for electron neutrinos. Note 
that the mean free paths in the central regions initially shrink as the star becomes denser 
and hotter, but by 15 s they begin to increase due to the decreasing temperature. It is also 
clear that the mu- and tau-neutrino mean free paths are more sensitive to the temperature 
than they are to the density, and their temperature dependence is greater than for electron 
neutrinos. The profiles for X Ufi are nearly mirror images of the temperature profiles for each 
time. Despite the fact that \ v /\ Ve > 1000 at early times in the core, there are so few mu- 
and tau-neutrinos relative to electron-neutrinos that the transport is dominated by v e s. 
However, by 15 s, X Vli /X Ue < 100 and the thermal evolution becomes dominated by mu- and 
tau-neutrinos. 
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5.2. Dependence Upon Initial Conditions 

A major uncertainty exists in the specifications of the initial entropy and Yl profiles. 
These are not accurately known because they depend on the details of the EOS and weak 
interaction physics employed both during the late evolution of the precollapse stars and the 
collapse of the core (Wilson & Mayle 1989, Swesty, Lattimer, & Myra 1994, Mezzacappa & 
Bruenn 1993). Ever since the work of Bethe et al. (1979), it has been realized that the inner 
core of the PNS should have canonical values of s ~ 1 and Yl ~ 0.38. Depending upon 
electron capture and beta decay rates and upon the symmetry energy, the entropy could be 
50% larger or smaller and Yl could vary by ±0.03. In addition, the s and Yl profiles in 
the shocked matter at enclosed baryon masses greater than 0.6 M (in the Mb = 1.6 M Q 
star) depend upon the strength of the outgoing shock and also, perhaps, upon convective 
instabilities. 

We performed several simulations of the Mb = 1.6 M Q GM3np star in order to study 
the sensitivity of the early evolution to the s and Yl profiles at t — 0. The spatial variation 
of various physical quantities for t — 0, 1, 5, and 15 seconds for models with lower (higher) 
initial central entropies are shown as dashed curves in Figure [TT] (Figure [12]), compared 
to the baseline model, which are displayed as solid lines. These variations do not produce 
dramatic effects, but the differences in thermodynamic quantities, although damped, persist 
to late times. Lowering the initial central entropy results in a higher value for since the 
central density is increased because of the lessening of thermal pressure. However, by 15 s, 
the values of \i v actually fall below that of the baseline simulation, because the evolution is 
more rapid due to the lower temperatures and relatively longer mean free paths. Exactly 
the opposite behavior occurs when the initial central entropy is increased. 



In Figure [13], we show the effects of increasing the lepton concentrations throughout 
the star; the central value of Yl is increased from 0.36 to 0.40. Since there are now more 
leptons to lose, the simulation results in more heating, which contributes toward longer 
evolution time scales. However, the effects are far smaller than in the case in which the core 
entropies were changed. 

We also considered variations of the entropy in the shocked mantle material; this case 



is shown as dotted lines in Figure |1^. Note that although the entropy is larger at all masses 
in the initial model, the temperatures are uniformly lower. This counterintuitive situation 
results from the additional expansion of the outer stellar layers: the density at all masses 
are also smaller in this case. Nevertheless, due to the rapid and extreme compression of the 
outer layers, these variations are completely swamped and there is virtually no difference 
between this and the baseline cases visible after a second. The same is true of simulations 
with variations in the lepton concentration in the shocked mantle material, which we do not 
display. Variations of entropy and lepton fraction in these outer regions manifest themselves 
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as short-term (a few tenths of a second) transients in the neutrino luminosities, which are 
discussed in §5.5. We can conclude that only variations of the initial conditions in the 
central regions of the PNS are relatively important in the long-term evolution. 

A useful summary of the effects of variations in initial conditions is contained in 



Figure [14], which displays the total neutrino luminosity as a function of time. Varying the 
central values of the entropy produce persistent differences in the luminosity. This is due 
to two factors: first, the initial binding energies of lower entropy models is larger, so the 
integrated luminosities must be smaller, and second, the evolution time scales of the lower 
entropy models must be more rapid, due to the lower temperatures and longer neutrino 
mean free paths. Similar behavior is produced by varying the central values of Y^. higher 
values of Y L reduce the initial binding energy and also lead to more entropy generation, 
higher temperatures, and longer diffusion times. Altering entropy and lepton profiles in the 
shocked mantle has little effect except in the first few tenths of a second. 

The observed variations in luminosity and average energy due to uncertainties in 
initial conditions are relatively small, but nevertheless will contribute to uncertainties in 
distinguishing EOS and compositional variations from an observed neutrino signal. 



5.3. EOS Dependence 

Equations of state differ in many respects at high densities. Most significant for the 
PNS are variations in compressibilities, symmetry energies, specific heats, and compositions. 
To study these effects, we consider three relativistic field-theoretical models with different 
high density properties. For two of these models, we also consider the possible presence of 
hyperons. These models do not necessarily encompass the entire range of possible behaviors. 
For example, we could have considered models with more extreme density dependences of 
the baryon effective masses, which are important in the determination of the specific heats, 
or the symmetry energy, which plays a role in determining the composition and the neutron 
star radius. However, the results of such variations may be gauged from the results we 
present here. 



The left panels of Figure [15| and Figure [16| compare the evolutions of a PNS of an 



Mb = 1.6 M Q star for the parameter sets GM3 and GM1 for nucleons-only matter for the 
1.6 M baryon mass star. GM1 is a stiffer EOS, and it has a lower value of the baryon 
effective mass M* . Therefore, the densities in the GM3np models are uniformly greater 
than in the GMlnp models. The GM3np temperatures are also greater; the higher densities 
more than compensate for the greater specific heat compared to GMlnp. As a result of 
the higher temperatures and densities, the neutrinos are more tightly trapped, there is 
more entropy generation, and the evolution times are longer in the GM3np models. When 
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hyperons are included, as shown in the right panels of Figure ^ and Figure [R], the GM1 
and GM3 models display the same relative trends. 

Comparing models with and without hyperons, i.e., the right and left panels of 
Figure [15] and Figure [16], the initial models (t = s) with hyperons show a small amount of 
softening. The central density of the hyperon stars becomes progressively larger than that 
of the purely nucleon stars, largely because of the appearance of thermal hyperons. The 
central densities of the hyperon stars remains larger even at late times when the remnant 
is cooling because hyperons are present. Little difference in the temperature evolution of 
these cases is apparent even after 20 s. Nevertheless, the evolution time scale of hyperon 
stars is slightly larger than for purely nucleon stars, due to the smaller mean free paths of 
hyperonic matter (see Figure |1(]). An interesting feature is the appearance of a trough in 
Y e and Y v near 0.75 M at t — 0.5 s. This feature is a short lived transient, because it does 
not appear in Figure [13| which shows Yl for t=l s. 

The major effects on the temporal evolution are more easily seen by examining 



Figure [L7|, which shows the time development of important physical variables at the star's 
center. It is clear that deleptonization continues during the subsequent cooling epoch. One 
should also note the plateau in the decay of Y v and \x v in the case of GM3np and GMlnp. 
This is due to the density sensitivity of the neutrino mean free paths. With hyperons, 
larger central densities and higher electron neutrino energies are reached. This increases 
the neutrino opacity, temporarily reducing the loss of neutrinos from the stellar core. 



Figure [18] shows the luminosities and the average neutrino energies as a function of time 
for the different EOS models, and for the baryon masses Mb = 1.6M and Mb = 1.8M . 
The luminosities are virtually identical until about 10 s, at which point the luminosity 
for models with a stiffer high density EOS decrease more rapidly. Contrasting the models 
with and without hyperons, we find a similar trend, with the hyperon models being able 
to sustain a higher luminosity at late times. The average energy of the emitted neutrinos 
are more sensitive to the EOS; once again a clear trend emerges in which models with a 
stiffer high-density EOS favor lower average energies. The average energy for models with 
hyperons is about 10-15% larger than those observed in matter containing only nucleons. 
These contrasting features are accentuated at late times and with increasing mass. 

We emphasize that for the mass (1.6 M ) and coupling constants chosen for this 



comparison relatively few hyperons are present at early times (see left panel of Figure [19]) . 
In fact, the maximum concentration of hyperons at t = is about 0.01%. However, as the 
density and temperature in the star is raised during deleptonization, the number and effect 
of hyperons becomes greater. Hyperons do not appear only at the star's center; the large 
temperatures in the mantle around 1 M result in a higher concentration of hyperons there 
than at the center from 0.25 s to 2 s. At late times, the central density of the npH star is 
nearly 20% larger than the np star and the concentration of hyperons at the star's center is 
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about 13%. For the mass (1.8 M ), a significantly larger number of hyperons are present as 
can be seen from the results shown in the right panel of Figure 19. The interesting feature 
to note is that at until about 15 s, the hyperons fractions are quite small, indicating that 
deleptonization in the inner core occurs very slowly. This is one of the reasons why we see 
very small changes in the observed luminosities between np and npH models at early times. 
In addition, feedbacks between the EOS and the opacities tend to reduce the differences. 



5.4. Mass Dependence 

We considered the effect of varying the initial mass of the PNS. For this study we 
employed the GM3 parameterization, and considered both np and npH matter. The actual 
baryon masses used in the simulations are listed in Table 3. The results displayed in 



Figures [0|, [2(J, and ^T] indicate that the most significant effect of increasing the initial mass 
is to increase the overall time scales of evolution, and to increase both the total and average 
energies released. The central heating also becomes more significant as one increases the 
remnant mass. Although the central entropy shows very little increase with increasing 
mass, the higher density of larger mass models translates into significant increases in central 
temperature. The combination of higher densities and temperatures lengthens the time 
constants associated with heating, cooling, and deleptonization. This naturally results in a 
longer duration of the plateau stage in the evolutions of Y v and \i v for np models. Note 
that these plateaus do not appear for npH matter for any mass. In spite of the larger time 
scales, the larger energy release, and increased temperatures larger mass models lead to 
larger luminosities at all times. 



5.5. Metastability 

The Mb = 1.8M baryon mass star containing hyperons is metastable. After more than 
100 s, this star will collapse to a black hole since the maximum baryon mass of the cold, 
catalyzed GM3npH EOS is 1.75 M . The possibility of metastability has been suggested in 
the context of SN 1987A, for which no neutron star remnant has been observed (Bethe & 
Brown 1995, Prakash, Cooke, & Lattimer 1995, Ellis, Lattimer, & Prakash 1996, Prakash et 
al. 1997). The lack of a neutron star in this object seems peculiar, because a neutron star 
rotating with the frequency of the Crab pulsar and with a magnetic field only a fraction 
of the Crab's field would be contributing noticeably to the light curve of the expanding 
gaseous remnant of SN 1987A. This is not observed. Although one cannot yet rule out the 
existence of a neutron star on the basis of Eddington-limit accretion, because the effective 
opacity may be due to iron instead of hydrogen and therefore be much greater, a black hole 
remnant seems a possibility. This hypothesis can be reconciled with the neutrino emission 
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observed from SN 1987A if the PNS was metastable because of the existence of hyperons or 
a kaon condensate in cold, catalyzed dense matter. 

For this reason, we further explored the case of a metastable PNS for the GM3 



parameterization. Figure |22| shows the time development of the central baryon density (top 
panel) and also the time to the collapse instability as a function of baryon mass (bottom 
panel). The larger the mass, the shorter the time to instability, since the PNS doesn't 
have to evolve in lepton number as much. Above 2.005 M , the metatstability disappears 
because the GM3npH initial model with the lepton and entropy profiles we chose is already 
unstable. Below about 1.73 M , there is no met ast ability, since this is the maximum mass 
of the cold, catalyzed npH star for GM3. These values are listed in Table 2, which contains 
corresponding, but larger, values for the case of GM1. The signature of neutrino emission 
from a metastable PNS should be identifiable and it is discussed in 55.7. 



Figure 123 shows the behavior of the quantity GMq/Rc as functions of both the radius 



and the enclosed baryon mass for a GM3npH model with Mb = 1.88 M Q . The time to 
instability for this model is slightly more than 45 s. Although the actual collapse to form 
a black hole is dynamic, and in any event could never be seen from the outside, one could 
estimate that the event horizon is likely to form where the redshift is largest. With time, 
it appears that this point moves outwards in the star, so that the event horizon can be 
expected to form relatively close to the star's physical surface and, therefore, outside the 
star's neutrinosphere. 



5.6. Convective Instabilities 

The possibility of convection in PNS has been discussed by a number of authors 
(Burrows & Fryxell 1993, Herant et al. 1994, Keil, Janka, & Miiller 1996, Mezzacappa et 
al. 1998). The decreasing entropy gradient in the star's exterior is obviously convectively 
unstable, and this gradient works its way into the interior on the deleptonization time scale. 
However, the complete determination of convective instability also requires knowledge of 
the lepton gradients. Explicitly, we consider the Ledoux criterion for convectively unstable 
regions: the condition is Cl(t) > 0, where 
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(65) 



We show in Figure the unstable regions of a PNS for the GM3np and GM3npH models 
for Mb = 1.6 M as shaded regions. Generally, the outermost few tenths of a solar mass is 
initially unstable. After a few seconds, the instability moves inwards, reaching the center 
after 10-15 seconds. We do not allow for convection in our simulations, so the result that 
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after a few seconds the exterior of the star becomes stable against convection is unreliable. 
In fact, the calculations of Keil, Janka, & Miiller (1996) show that the convection proceeds 
similarly to what our models show, but that the exterior remains convectively unstable 
throughout the evolution. Their time scales are, however, much shorter in comparison 
to our findings. The influence of convection upon the luminosities and average neutrino 
energies has been shown by prior authors (Burrows & Lattimer 1988) to be relatively 
significant, and we will consider this effect in the next generation of our calculations. 



5.7. Signals in Terrestrial Detectors 

Knowledge of the neutrino luminosities and energies permit the signal from a supernova 
to be calculated. To be specific, we will present results in a way that makes comparison to 
the signal observed from SN 1987A (Bionta et al. 1987, Hirata et al. 1987) transparent. We 
assume two fiducial detectors, KII and 1MB, whose masses, thresholds and efficiencies are 
taken to be the same as in Lattimer & Yahil (1989). We focus on the antineutrino signal, 
assuming that the antineutrino luminosity L v is 1/6 of the total. 

We assume a pure water detector for which the cross section of antineutrinos absorbed 
onto protons is cr(E u ) = oo-E 2 where E v is the antielectron neutrino energy in MeV and 
o"o = 9.3 ■ 10~ 44 cm 2 . There are n p = 6.7 ■ 10 28 free protons per kiloton of water. We take 
the distance to SN 1987A to be D = 50 kpc. Assuming that the neutrinos leaving the 
protoneutron star have a Fermi-Dirac spectrum with zero chemical potential, the count rate 
is 

f = (f )' ^ M £ E ^ E - • «*» 

where E t h is the detector threshold, W is the detector efficiency and M. is the detector mass 



in kilotons. Also, = e~^R v and T£° = e?'T v , where e* 8 = y 1 - 2GM G /R u c 2 and M G 
and R v are the gravitational mass of the neutron star and the radius of the neutrinosphere, 
respectively. Since the luminosity is calculated more precisely than the neutrinospheric 
radius and temperature, we rewrite this, substituting the antineutrino luminosity, as 

dN a n G,(E th ,e^T u ) 

dt ~ AnD^ F 3 (0) v v > [ } 

where -F 3 (0) = 77r 4 /120 is an ordinary Fermi integral and Gi(E t h,T) denotes a modified, 
truncated Fermi integral: 

POO 

Gi(E th ,T)= dzz*W(zT)(l + e z )- 1 ■ (68) 



th 



In Figure |25| we display the count rate and integrated counts for KII and 1MB for 
the M B = 1.6 and 1.8 M models for the various EOS cases studied. The results can be 
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easily interpreted in terms of the average energies illustrated in Figure [18|. While the total 
luminosities remain essentially identical for all cases of a given mass up to about 20 seconds, 
the average energies show a systematic variation. The softer the underlying EOS, the larger 
the average energy, and this directly translates into higher count rates. The general trends, 
which are explainable by the more extreme temperature sensitivity of the 1MB detector, 
are that more counts accumulate in KII than in 1MB and the time over which a signal is 
observed is also larger for KII than for 1MB. 

The total variation among the EOS models in the accumulated counts after 40 seconds 
is about 10% for a given mass. The mean ratio of total counts between the 1.8 M and 1.6 
M cases is about 1.5. Therefore, it seems that for a Galactic supernova observed in the 
near future, in which the expected number of counts will be substantially larger because of 
its proximity and larger detectors, that the protoneutron star mass could be measured to 
better than 0.05 M from the total number of counts observed. It also seems that the more 
subtle differences between EOS models will be observable only if the accumulations after 
about 20 seconds are large enough. 



5.8. Comparison with Earlier Work 

We turn now to a comparison of our results with the earlier works of BL and KJ 
on the early evolution of protoneutron stars containing only nucleonic matter. There are 
important differences arising due to the different input physics. In these earlier works the 
neutrino opacities were not calculated consistently with the dense matter EOS. Thus, it is 
possible to compare only some aspects of previous calculations with our own. These aspects 
include the luminosity evolution and the average neutrino energies, as well as the evolution 
of temperature and neutrino number in the interior. 

BL and KJ both found that, for 1.4-1.6 M Q star, the deleptonization time, the time 
needed for the electron neutrino number Y Ve to vanish in the star's center, was about 10 
s (see Figure 8 in BL & Figure 5 in KJ). In addition, they both found that the central 
temperature reached a maximum at about the time at which the neutrino number vanished. 
In contrast, we generally find that these two time scales are not similar: while the time 
required for the central temperature to reach a maximum is of order 10 seconds, the 
deleptonization time is much longer, at least 30 seconds. We attribute the more rapid 
deleptonization in earlier works to inadequate treatments of the temperature dependence of 
the direct Urea process, the nuclear symmetry energy and the matter's specific heat. In our 
models, the larger electron neutrino cross sections at high temperature and low neutrino 
numbers stem the flow of electron neutrino number. 

Comparing the calculations of BL in more detail, we find that the total luminosity is 
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larger during the deleptonization era. After 20 seconds in BL, nearly 1/3 of the binding 
energy has yet to be radiated (see Figure 19 in BL), while in our calculations at the same 
time less than 10% remains. Burrows, Klein, & Gandhi (1992) have argued that the 
neutrino luminosity evolves as a power law at late times (t larger than a few seconds), 
while it decays exponentially in ours. As Prakash et al. (1995) demonstrated, an effective 
E 2 behavior of the neutrino cross sections should result in an exponential decrease of 
temperature and luminosity for the star. 

Another apparent difference between the calculations of BL and ours concerns the 
average emitted neutrino energy: in our calculation this peaks at about 5 s, while in BL 
it steadily decreases with time (see Figure 20 in BL). This difference could be important 
for the questions of neutrino-assisted shocks and for nucleosynthesis. However, it should be 
borne in mind that we estimated the average energy by the relation < E u >= 3T„, but this 
neglects possible effects of neutrino degeneracy and the fact that the emergent spectrum 
may not be Fermi-Dirac. It will be necessary to employ more accurate transport in the 
vicinity of the neutrinosphere before firm conclusions can be drawn. 

We see more substantial differences from the calculations of KJ. During the first few 
seconds of the evolution, KJ find that the central density decreases (see Figure 10 in KJ). 
KJ attribute this decrease to the heating which occurs during deleptonization. Indeed, KJ 
find a large increase in central temperature during the first five seconds (see Figure 11 in 
KJ) which is much more rapid than in our calculations. The effects of rapid deleptonization 
and rapid temperature increase could be explained by a charged- current mean free path 
that is several times larger than in our calculation. On the other hand, the cooling time 
scales, which are controlled by neutral current processes, are rather similar to what we have 
found. This probably also accounts for the remarkable similarities between our and KJ's 
total neutrino luminosities, despite very different delptonization times. 

A decrease in central density during the deleptonization epoch could happen if the 
thermal contributions to the pressure at 3-5 times the nuclear density were a significant 
part of the total. However, this would imply that the specific heat of the matter is radically 
different than what we find. 

We can also compare our protoneutron star evolutions to KJ for cases in which 
hyperons appear. In contrast to their calculations, we find only minor differences in the 
evolution of nucleons-only and hyperon-bearing models at early times. In fact, for most 
masses studied, the evolutions of luminosity in our calculations for the two cases are nearly 
identical for about 20 seconds. KJ find that large differences in the luminosity occur 
between np and npH cases as early as 2-4 seconds (see Figure 9 in KJ) in the evolution. 
They found a time to collapse in the metastable case of less than 20 seconds when the 
stellar mass was 0.02 M larger than the critical mass (see Table 2 in KJ). For the same 
mass difference, we find a time to collapse of more than 100 seconds. These differences no 
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doubt stem from the same effects that led to the large differences apparent between KJ's 
and our nucleon models. Our relatively large deleptonization times prevent hyperons from 
flooding the interior until rather late times. This significantly increases the time spent in 
the metastable phase when the mass is only marginally above the maximium mass for cold, 
deleptonized matter (see Figure p2j). 

6. Summary and Outlook 

In this paper, we presented calculations of protoneutron star evolutions which used 
neutrino opacities that were consistently determined from the underlying dense matter 
EOS. We explored the sensitivity of the results to the initial stellar model, the total mass, 
the underlying EOS, and the addition of hyperons. We find that the differences in predicted 
luminosities and emitted neutrino energies do not depend much upon the details of the 
initial models or the underlying high-density EOS for early times (t < 20 s), provided that 
opacities are calculated consistently with the EOS. The same holds true for models which 
allow for the presence of hyperons, except when the initial mass is significantly larger 
than the maximum mass for cold, catalyzed matter. Increasing the total mass tends to 
increase luminosities and average energies at all times. For times larger than about 10 
seconds, and prior to the occurrence of neutrino transparency, the neutrino luminosities 
decay exponentially with a time constant that is sensitive to the high-density properties 
of matter. We also find the average emitted neutrino energy increases during the first 5 
seconds of evolution, and then decreases nearly linearly with time. Metastable stars, those 
with hyperons which are unstable to collapse upon deleptonization, have relatively long 
evolution times, which increase the nearer the mass is to the maximum mass supported by a 
cold, deleptonized star. 

There are several interesting and possibly important physical effects that remain to be 
investigated. 

(1) In terms of input neutrino opacities, this work provides a baseline calculation in terms 
of which further improvements (Reddy et al. 1997, Prakash & Reddy 1997, Burrows & 
Sawyer 1998a;b, Hannestad & Raffelt 1998) such as the effects of many-body correlations, 
axial charge renormalization and multi-pair excitations may be assesed. 

(2) While the diffusion approximation is valid for a significant portion of the star at early 
times, neutrino transport in the semi-transparent regions is better handled in a multigroup 
transport scheme. Better transport in the semi-transparent regions will play an important 
role in determining the average neutrino energies for the different neutrino species and their 
energy spectra. 

(3) An adequate treatment of convection coupled with neutrino transport appears to be 
necessary based on the large regions we identified as being convectively unstable. 
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(4) Other softening components in dense matter might produce effects dissimilar to those 
we find when considering hyperons. For example, the appearance of a Bose condensate 
differs from the appearance of hyperons because the condensate does not contribute to the 
matter's specific heat. In addition, the presence of a condensate could produce distinctive 
modifications to the neutrino opacities. 

Results of calculations incorporating these improvements will be reported subsequently. 
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TABLE 1 



NUCLEON-MESON COUPLING CONSTANTS 



Model 


K 


O'sym 


M*/M 


9 a /m a 




9 P /m p 


b 


c 


GM1 


300 


32.5 


0.70 


3.434 


2.674 


2.100 


0.002950 


-0.00107 


GM3 


240 


32.5 


0.78 


3.151 


2.195 


2.189 


0.008659 


-0.002421 


GM4 


240 


27.5 


0.78 


3.151 


2.195 


1.862 


0.008659 


-0.002421 



NOTE.- Coupling constants for the GM (Glendenning & Moszkowski 1991) models which 
are calibrated at the nuclear equilibrium density no = 0.153 fm~ 3 to give an energy 
per particle of symmetric nuclear matter of —16.3 MeV. The compression modulus K Q , 
symmetry energy a sym and Dirac effective mass M* for the different models are also listed. 

TABLE 2 



MAXIMUM GRAVITATIONAL (BARYONIC) MASSES 



Model 


Y L = 0.4 (s = 1) 


Y v = (s = 2) 


Y v = (s = 0) 


GMlnp 
GM3np 
GMlnpH 
GM3npH 


2.278(2.609) 
1.946(2.181) 
2.024(2.259) 
1.768(1.944) 


2.367(2.748) 
2.044(2.328) 
1.793(1.991) 
1.573(1.726) 


2.346(2.805) 
2.005(2.345) 
1.776(2.020) 
1.544(1.732) 



NOTE.- Symbols are Yl: the total lepton fraction, Y Ve : electron neutrino fraction, and s: 
the total entropy per baryon (in units of Boltzmann's constant). Model designations are as 
in Table 1 with np referring to nucleons-only matter and npH referring to matter including 
hyperons. 

TABLE 3 



BARYONIC MASSES OF PNS EVOLUTION MODELS 



Mass Label 


Actual Mass (M ) 


Mass Label 


Actual Mass (M ) 


1.08 


1.0852 


1.84 


1.8434 


1.4 


1.4038 


1.88 


1.8843 


1.6 


1.6125 


1.90 


1.9038 


1.8 


1.7981 


1.94 


1.9413 


1.81 


1.8162 


1.98 


1.9863 
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FIGURE CAPTIONS 



FIG. 1.- Relative concentrations in lepton-poor (left panels) and lepton-rich (right panels) 
matter at finite entropy as a function of baryon density % (no is the nuclear equilibrium 
density). Top panels show results in nucleons-only matter. The bottom panels are for 
matter including hyperons. 



FIG. 2.- Total pressure including contributions from leptons in lepton-poor (left panels) 
and lepton-rich (right panels) matter at finite entropy as a function of baryon density % 
(no is the nuclear equilibrium density). The symbol np refers to matter with nucleons-only 
and npH to matter including hyperons. 



FIG. 3.- Temperatures attained in lepton-poor (left panels) and lepton-rich (right panels) 
matter at finite entropy as a function of baryon density % (no is the nuclear equilibrium 
density). The symbol np refers to matter with nucleons-only and npH to matter including 
hyperons. 



FIG. 4.- Lepton chemical potentials in lepton-poor (left panel) and lepton-rich (right panel) 
matter at finite entropy as a function of baryon density % (no is the nuclear equilibrium 
density). The symbol np refers to matter with nucleons-only and npH to matter including 
hyperons. 



FIG. 5.- Binding energy as a function of baryon mass. The symbol np refers to matter with 
nucleons-only and npH to matter including hyperons at T=0 and Y v =0. 



FIG. 6.- Neutrino mean free paths in matter with nucleons only (left panels). Right panels 
show ratios of mean free paths in matter without and with hyperons. Abscissa is baryon 
density n B (n is the nuclear equilibrium density). Top panels show scattering mean free 
paths common to all neutrino species. The bottom panels show results for electron neutrino 
mean free paths where absorption reactions are included. The neutrino content is labelled 
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in the different panels. 



FIG. 7.- Inverse neutrino mean free paths due to charged current reactions versus 
temperature. 



Fig. 8.- Comparison of scattering mean free paths in neutrino poor matter at fixed entropy 
for different EOSs in matter containing nucleons and also hyperons. 



FIG. 9.- Temporal variation of thermodynamic quantites (s: the entropy per baryon, 
T: the temperature, the baryon density Y u : the net electron neutrino fraction, \i v : 
electron neutrino chemical potential, and Y e : the net electron concentration) inside a star 
for the baseline GM3 model for nucleons-only matter. Labels indicate time in seconds. The 
abscissa shows the enclosed baryonic mass. 



FIG. 10.- Typical electron-neutrino and muon-neutrino mean free paths. For the electron 
neutrino the average energy < E v >= \L Ve + ttT, while for the muon and tau neutrinos 
< E v >= 7rT. Top panels show temporal variations for the baseline GM3 model in 
nucleons-only matter and the bottom panels show results in matter including hyperons. 



FIG. 11- Effects of changing initial entropy profiles. Results of the baseline GM3 model 
(solid curves) with nucleons only with an inital s ~ 1 in the core are compared to those 
(dashed curves) of a model with reduced initial entropy in the core (s ~ 0.5). Symbols as 
in Figure [|. 



FIG. 12.- Effects of changing initial entropy profiles. Results of the baseline GM3 model 
(solid curves) with nucleons only with an inital s ~ 1 in the core are compared to those 
(dashed curves) of a model with an enhanced initial entropy in the core (s ~ 2). The dotted 
curves show effects of enhancing the entropy in the surface regions. Symbols as in Figure |9|. 
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FIG. 13.- Effects of changing initial lepton profiles. Results of the baseline GM3 model 
(solid curves) with nucleons only with an inital Y L = 0.38 in the core are compared to those 
(dashed curves) of a model with an enhanced initial lepton content in the core (Yj, = 0.4). 
Symbols as in Figure [5| 



FIG. 14.- Comparison of total neutrino energy luminosity versus time between the baseline 
model and different initial models. 



FIG. 15.- Temporal variation of the entropy, temperature and baryon density in matter 
containing nucleons only (left panels) and in matter that also contains hyperons (right 
panels). Differences arising from the underlying EOS are apparent by comparing the solid 
and dashed curves. The abscissa shows the enclosed baryonic mass in a 1.6M star. 



FIG. 16.- Temporal variation of the net electron neutrino fraction, electron neutrino 
chemical potential, and the net electron concentration in matter containing nucleons only 
(left panels) and in matter that also contains hyperons (right panels). Differences arising 
from the underlying EOS are apparent by comparing the solid and dashed curves. The 
abscissa shows the enclosed baryonic mass in a 1.6M star. 



FIG. 17.- Effects of changing the underlying EOS (GM1 versus GM3) on the core values of 
the thermodynamic quantities versus time in matter containing nucleons only (thin solid 
and dashed curves) and in matter that also contains hyperons (thick solid and dashed 
curves). Results are for a M B = 1.6M star. 



FIG. 18.- Comparison of the mean neutrino energy (top panels) and the total neutrino 
energy luminosity (bottom panels) versus time for the different EOS models. The left 
panels show results for a baryon mass Mb = 1.6M Q and the right panels for a baryon mass 
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M B = 1.8M ( 



FIG. 19.- Temporal evolution of strangeness per baryon for the indicated baryon masses. 



FIG. 20.- Mass dependence of the core values of the thermodynamic quantities versus time 
in matter containing nucleons only (solid curves) and in matter that also contains hyperons 
(dashed curves). 



FIG. 21.- Mass dependence of the average energy of the emitted neutrinos (top panels) and 
the total neutrino energy luminosity (bottom panels) in models with and without hyperons. 



FIG. 22.- Top panel: Evolution of the central baryon number density for different baryonic 
mass stars containing hyperons (model GM3npH) which are metastable. Bottom panel: 
Time required by stars shown in the top panel to reach the unstable configuration. 



FIG. 23.- Evolution of (GM G /Rc 2 ) in a M B = 1.88M star containing hyperons (model 
GM3npH). Inset labels show time in seconds. 



FIG. 24.- Evolution of convectively unstable regions in a M B = 1.6M star. Top panel 
shows results for a star containing nucleons only and the bottom panel is for a star that 
also contains hyperons. 



FIG. 25.- Neutrino signals in terrestrial detectors KII and 1MB, and for a supernova at 50 
kpc. Top panels show total number of counts, and the bottom panels show the count rate. 
Left panels are for M B =1.6 M , and right panels show results for M B =1.8 M . 
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Fig. 6.- 
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Fig. 7.- 
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Fig. 9.- 
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Fig. 12.- 




Fig. 13.- 
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Fig. 14.- 
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Fig. 17.- 
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Fig. 18.- 
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Fig. 19.- 
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Fig. 25.- 



